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ABSTRACT 

We investigate the process of formation and subsequent evolution of promi- 
nence plasma in a filament channel and its overlying arcade. We construct a 
three-dimensional time-dependent model of a filament- channel prominence suit- 
able to be compared with observations. We combine this magnetic held structure 
with one-dimensional independent simulations of many flux tubes. The mag- 
netic structure is a three-dimensional sheared double arcade, and the thermal 
non-equilibrium process governs the plasma evolution. We have found that the 
condensations in the corona can be divided into two populations: threads and 
blobs. Threads are massive condensations that linger in the held line dips. Blobs 
are ubiquitous small condensations that are produced throughout the filament 
and overlying arcade magnetic structure, and rapidly fall to the chromosphere. 
The total prominence mass is in agreement with observations. The threads are 
the principal contributors to the total mass, whereas the blob contribution is 
small. The motion of the threads is basically horizontal, while blobs move in all 
directions along the held. The peak velocities for both populations are compara- 
ble, but there is a weak tendency for the velocity to increase with the inclination, 
and the blobs with motion near vertical have the largest values of the velocity. 
We have generated synthetic images of the whole structure in an Ha proxy and 
in two EUV channels of the AIA instrument aboard SDO. These images show the 
plasma at cool, warm and hot temperatures. The theoretical differential emis- 
sion measure of our system agrees very well with observations in the temperature 
range log T = 4.6 — 5.7. We conclude that the sheared-arcade magnetic structure 
and plasma dynamics fit well the abundant observational evidence. 


1 CRESST and Space Weather Laboratory NASA/GSFC, Greenbelt, MD 20771, USA 

2 NASA/GSFC, Greenbelt, MD 20771, USA 

3 Naval Research Laboratory, Washington, DC 20375, USA 



1. Introduction 


Quiescent solar prominences have been studied for a long time, and significant progress 
in determining their nature has been made. We can say with certainty that they are com- 
posed of cold plasma, and it is generally accepted that the magnetic field supports these 
dense structures in the tenuous corona. However, there are still important unknowns about 
their structure, formation and evolution. High-resolution observations from space and from 
the ground - SOHO, TRACE, SVST, VAULT, Hinode, and SDO - show a complex fine 
structure in space and time. This fine structure consists of moving blobs and bundles of 
horizontal threads with counterstreaming motions (Engvold 1976; Zirker et al. 1998; Vourl- 
idas & Korendyke 2003; Lin et al. 2003; Lin et al. 2005; Berger et al. 2008). The quiescent 
prominences are immersed in structures even larger than themselves. Prominences form 
above a polarity inversion line (Babcock & Babcock 1955; Smith & Ramsey 1967; McIn- 
tosh 1972) in filament channels (Martin 1998; Gaizauskas et al. 2001), which are regions 
of sheared magnetic field (Hagyard et al. 1984; Moore et al. 1987; Venkatakrishnan et al. 
1989). Their observational signature in the chromosphere is fibrils forming a path aligned 
with the polarity inversion line (Smith 1968; Praia 1971; Foukal 1971; Rompolt & Bogdan 
1986; Tandberg-Hanssen 1995). Magnetic loop arcades overlie filament channels, spanning 
high above them and rooted at both sides of the channel (Tandberg-Hanssen 1995; Martin 
1998). Typically these loops are at the base of a helmet streamer (Engvold 1989). The over- 
lying arcade is visible in white light, in soft X-rays, and in UV lines. If the filament is viewed 
from one end, a dark coronal cavity frequently can be observed between the prominence and 
the overlying arcade (Tandberg-Hanssen 1995; Martin 1998). Direct measurements of the 
magnetic fields have shed light on the magnetic structure of quiescent prominences (Leroy 
1989; Casini et al. 2003; Kuckein et al. 2009). However, these measurements are focused on 
the cool prominence plasma alone and do not establish definitively the global structure of 
the magnetic field in the coronal portion of the filament channel. 

There are basically two models of the magnetic structure of filament channels: the 
sheared arcade and the flux rope (see review by Mackay et al. 2010). Antiochos et al. 
(1994), DeVore & Antiochos (2000), and Aulanier et al. (2002) showed that the sheared 
arcade model exhibits most of the observed properties of the large-scale filament channel 
structure. The authors consider that the magnetic structure is a natural consequence of 
observed photospheric shear motions along polarity inversion lines (see, e.g., Welsch et al. 
2004). The generated magnetic field structure has dipped and flattened field lines consistent 
with most prominence observations (Martin 1998). Many quiescent filaments are observed to 
grow by linkage between smaller segments (Malherbe 1989; Gaizauskas et al. 2001; Schmieder 
et al. 2004a; van Ballegooijen 2004). DeVore et al. (2005) and Aulanier et al. (2006) studied 
the interaction of two neighboring sheared bipoles with different combinations of polarity 
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and chirality. They found that merging two adjacent arcades with aligned axial field along 
their spines formed a long, stable, strongly sheared prominence. Assuming that the cool 
plasma resides in magnetic clips, the resulting structure reproduces the global appearance of 
a typical low-latitude intermediate quiescent filament (Mackay et al. 2010). 

It is generally believed that the origin of the prominence plasma is the chromosphere 
(Pikel’ner 1971; Saito & Tandberg-Hanssen 1973; Zirker et al. 1994), but the mechanism that 
puts the chromospheric material into the corona is unknown. The thermal non-equilibrium 
process has been studied extensively (Antiochos & Klimchuk 1991; Antiochos et al. 1999, 
2000; Karpen et al. 2001, 2003, 2005, 2006; Karpen & Antiochos 2008), and yields quan- 
titative predictions of prominence plasma directly comparable with observations. In this 
model a localized heating in the flux tube footpoint produces an evaporation of the chro- 
mospheric plasma, which condenses in the coronal part of the tube. The localized heating 
near the chromosphere is consistent with several coronal-loop studies (Aschwanden et al. 
2001; Schmieder et al. 2004b). To produce this instability the heating must be concentrated 
in a region of approximately one eighth of the tube length or less (Serio et al. 1981; Mok 
et al. 1990; Antiochos & Klimchuk 1991; Dahlburg et al. 1998). In our earlier investiga- 
tions of the thermal non-equilibrium process, we found that large condensations form along 
long, low-lying magnetic field lines. In most cases the plasma undergoes repeated dynamic 
cycles of plasma condensation, streaming along the tube, and falling to the chromosphere. 
Small, short-lived knots are produced in coronal loops, possibly explaining the coronal rain 
(Schrijver 2001; Miiller et al. 2003, 2004, 2005). In deeply clipped flux tubes condensations 
form readily. The subsequent evolution depends on the tube slope. For slopes steeper than 
a critical value determined by the heating asymmetry and gravitational scale height, the 
condensations fall to the bottom of the dip and remain there accreting mass continuously. 

In the present work, we built a comprehensive model of a multi-threaded prominence 
combining the results of the three-dimensional sheared double-arcade model of DeVore et al. 
(2005) for the magnetic structure with the thermal non-equilibrium process for the plasma 
evolution in many independent one- dimensional simulations. This three-dimensional model 
allows us to study the basic features of a multi-stranded prominence and the process of 
formation and evolution of its plasma. This paper is arranged as follows. In §2 we describe 
the theoretical model used in the current work. The detailed results of the investigation 
are presented in §3. The global appearance of the prominence in different spectral channels 
and the theoretical differential emission measure are shown in §4. Finally, in §5 the results 
of this investigation are summarized and conclusions are drawn. In particular, we conclude 
that our three-dimensional model fits well the observational evidence. 



2. Theoretical model 


As in our earlier studies, we assume that the magnetic structure is static on the timescale 
of the plasma evolution (see e.g., Martin 1998) and that the gas pressure of the plasma is 
negligible in contrast with the magnetic pressure. Additionally, in the solar corona the 
conductive heat flux is mainly constrained to follow the field lines. With all these con- 
siderations each flux tube is a thermally isolated system where the plasma motions are 
along the magnetic held. The equations describing the plasma evolution are reduced to an 
one-dimensional set of hydrodynamic equations for mass, momentum and energy conserva- 
tion (see, e.g., Karpen et al. 2005) including the Raymond-Klimchuk radiative loss function 
(Klimchuk & Cargill 2001), Spitzer (1962) thermal conductivity, and a prescribed volumet- 
ric heating. Hence the time-dependent plasma evolution along each static held line can be 
determined by solving numerically those equations. 

We use the magnetic held structure given by the sheared double-arcade model of DeVore 
et al. (2005) in Cartesian geometry. This magnetic held configuration is the result of shearing 
two identical dipoles with the same polarity, the same chirality, and with their axes aligned. 
The shear is produced by imposing differential footpoint motions over the photospheric 
plasma near the polarity inversion line (PIL). As a consequence, an inner bundle of sheared 
held lines is produced (see Fig. 1). These are long, low-lying held lines nearly parallel to the 
PIL, and many of them have dips. An outer bundle of non-sheared held lines comprises the 
overlying arcade with arched held lines perpendicular (or slightly skewed from perpendicular) 
to the PIL. In this model, the filament channel (FC) magnetic configuration is a natural 
consequence of the magnetic shear in a 3D structure. However, the exact origin of the 
magnetic shear is not important in our model. The shear could be the result of the partial 
emergence of hux tubes from below, or direct subsurface motions near the PIL. 

We performed the numerical simulations of the one-dimensional hydrodynamic equations 
with our code ARGOS (Adaptively Refined GOdunov Solver) with adaptive mesh refinement 
(see Antiochos et al. 1999, for details). We have selected a discrete set of 125 non-uniform 
hux tubes of the 3D sheared double- arcade magnetic structure shown in Figure 1. We have 
rescaled the computational domain of DeVore et al. (2005) to a box of 800 Mm parallel 
to the PIL, 200 Mm perpendicular to this PIL in the base plane, and 200 Mm in the 
vertical direction. The footpoints of the held lines were selected by imposing two uniform 
rectangular grids in the photosphere. The hrst grid is denser and covers the sheared region 
of the photosphere, being 40 Mm wide centered at the PIL and 600 Mm along the PIL. 
The second grid covers the outer region with no shear, and it is located only on one side 
of the PIL due to the symmetry of the magnetic held with respect to the PIL. This grid 
extends perpendicular to the PIL from the end of the hrst grid to 55 Mm away. Along the 
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Fig. 1. — Magnetic field structure of the sheared double arcade in Cartesian geometry. White 
lines are the magnetic field lines we have used in this work. In (a) the magnetic structure is 
viewed from above. At the bottom of the corona the normal magnetic field is plotted (color 
scale). The two dipoles are clearly seen in this view. In (b) the side-view orientation is 
shown. Finally, in (c) the end-view orientation is plotted. 

PIL both grids have the same extension. The geometry of the 3D flux tubes is incorporated 
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Fig. 2. — Flux tubes can be roughly classified into three groups according to their shapes. 
Here we show two examples of each group. The horizontal d axis connects the two tube feet. 
Flux tubes with a deep dip (solid lines) are sheared lines, and form the filament channel 
core. Shallow-dip tubes (dotted lines) are also sheared lines, and cross the polarity inversion 
line above the deep-dip lines. Finally, the loop-like tubes (dot-dashed lines) have no dip, 
and form the overlying arcade structure. 

into our simulations. The field line shape determines the projection of the gravity along the 
tube. In addition the non-uniform cross-sectional area A is incorporated using conservation 
of magnetic flux along the tube. The position s along the tube is described by the distance 
from one footpoint at the base of the corona. With this definition the positions of the two 
footpoints at the base of the corona are s = 0 and s = L, where L is the length of the coronal 
part of the tube. Setting the area at s = 0 to be H 0 , we can determine the area along the 
entire tube using A(s)B(s) = A 0 B 0 where B 0 is the magnetic field strength of the sheared 
double arcade at s — 0. 

The geometry of the flux tubes can be classified roughly into three groups in the whole 
magnetic structure: tubes with deep dips, tubes with shallow dips, and tubes with no dips 
(loop-like tubes). In Figure 2 two examples of each group are plotted. The first group has 
deep dips; these are sheared, long, and low-lying lines with inverse polarity at their central 
parts (tubes 1 and 2 in Fig. 2). The magnetic field and gravity jointly provide a stable 
equilibrium region in the dipped portion of the tube because the projected gravity forces the 
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plasma to move to the center of the clip. The second kind of tube has a plateau shape with 
a central part almost horizontal, as we see in the examples 3 and 4 of Figure 2, providing 
a region in the tube where the projected gravity is small or zero. These two groups of held 
lines belong to the sheared magnetic held, and they form the FC structure (Fig. 1). The 
third kind of tube is in the unsheared or weakly sheared magnetic held arcade overlying the 
FC. These tubes are loop-shaped as exemplihed by tubes 5 and 6 in Figure 2. Most of the 
loops cross the PIL perpendicularly, but the short tubes near the FC structure are skewed 
from the perpendicular direction. Most loop- like tubes have a large projected gravity except 
in the tube apices. There is a small group of loop-like tubes at both ends of the FC with 
large inclinations (see Figs, la and lb), with small projected gravity along the tube. Tube 
5 is a short loop that wraps the inner core of the FC structure and is slightly skewed from 
perpendicular direction to the PIL. Tube 6 is a very long loop reaching very high altitudes 
of around 190 Mm. Most of the selected 125 tubes belong to the FC, and the others form 
the overlying arcade. 

In Figure 3 the histogram of the lengths of the flux tubes is plotted. The distribution 
has a Gaussian shape except for the contribution of very long tubes. The tube lengths are 
distributed between 84 Mm and 479 Mm with a mean value of 218 Mm. We compute also 
the contribution of the three kinds of flux tubes to the global distribution of tubes. We 
see that the right wing of the whole distribution is given mainly by the tubes of the first 
kind with deep dips; the main contribution to the left wing of the distribution comes from 
both shallow-dip and loop held lines. Thus, on average the tubes with deep dips are longer 
than the tubes with shallow dips. This is because lines in the first group are reconnected 
in the process of formation of the FC by shear of the footpoints (DeVore et al. 2005). We 
have selected only a small number of tubes of the overlying arcade, as can be seen in the 
histogram. This is because we are interested mainly in the tubes of the FC core, but we 
also included a small sample of other tubes to take into account the evolution of the plasma 
outside the FC. The area expansion, A(s)/k4 0 , of the tubes of the first and second kinds in 
the corona is similar, with a mean factor of approximately 3 in the central part of the tubes. 
In contrast, the area expansion of tubes of the overlying arcade has an average value of 30 
in the loop apices. 

As in our previous studies, the volumetric heating has two components: a spatially 
uniform background heating, Hq, plus a localized component at each footpoint in the chro- 
mosphere, H\ and H 2l that falls off exponentially in the corona. Thus, 

H(s) =H 0 + Hi e~ s ' x + H 2 e (*- L )/\ (1) 

where A is the heating deposition scale that we have set to A = 10 Mm, and the background 
heating is set to H 0 = 1.5 x 10~ 5 erg cm -3 s _1 . The tube feet are heated asymmetrically, so H 1 




Fig. 3. — Histogram of the lengths of the selected flux tubes of the magnetic structure. The 
full set of 125 (black dashed line) has a distribution similar to a Gaussian. The minimum 
length is 84 Mm and the maximum is 479 Mm, with a mean value of 218 Mm. Tubes of 
the first kind (blue line) with deep dips have a Gaussian distribution with a mean value 
of 255 Mm. The second kind (green line), tubes with shallow dips, also have a Gaussian 
distribution but centered at 156 Mm. The third kind (red line), the overlying arcade loops, 
are tubes with a mean value of 192 Mm. 

and H 2 are different. We have set the larger value of the volumetric heating as max(H \ , H 2 ) = 
10~ 2 erg cm -3 s -1 , and the heating asymmetry as min(i7i, H 2 )/max(Hi, Ho) = 0.75. In 
our 125 independent simulations we have distributed randomly the position of the larger 
localized heating, so some tubes are more heated at s = 0 and others at s = L. A 60 Mm 
chromospheric region has been added to each flux tube footpoint in order to have a mass 
reservoir and sink for the coronal plasma evolution. In all flux tubes, an equilibrium situation 
was established first with only background heating turned on during the first 10 5 s. This 
initial uniform heating populates the flux tubes with density, temperature, and pressure 
dictated by the scaling laws of Rosner et al. (1978). After this 10 5 s, the localized heating 
was ramped up over 1000 s and sustained at that level thereafter. The simulation time is 
10 5 s from the onset of the localized heating (t = 0). Additional details of the numerical 
technique can be found in Antiochos et ah (1999) and Karpen et ah (2005). 
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3. Results 

In this work, we consider the full system of tubes and the features of the whole multi- 
stranded prominence. The evolution of the plasma on each flux tube depends on the tube 
geometry and the asymmetry of the heating function (Karpen et al. 2003). As we have seen 
in §2 the geometry of the field lines can be classified roughly into three groups. In Figure 
2, two representative tubes of each group have been plotted. In order to understand the 
process of formation of the whole prominence, we first study the evolution of the plasma on 
the selected set of six field lines. In Figure 4, we have plotted the mass of the condensations 


(a) 



(b) 



(c) 



time (hrs) 


Fig. 4. — Temporal evolution of the condensation masses per unit area for the six represen- 
tative tubes, where the masses are divided by A 0 . (a) Condensation mass of the two tubes 
with deep dips, (b) Condensations formed in tubes with shallow dips, (c) Condensations 
in the two overlying arcade tubes. 

as a function of time for the six simulation examples. Figure 4a shows the evolution of 
cool condensations in tubes 1 and 2 with deep dips. The condensation in tube 1 forms at 
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t — 3 hrs after the onset of the localized heating, after which the condensation mass grows 
linearly at a rate of 6.6 x 10 -4 g cm -2 hr -1 for the rest of the simulation. The condensation 
forms at the dipped portion of the held line and remains there while it accretes mass. The 
condensation of tube 2 shows completely different behavior, forming close to a footpoint and 
falling rapidly to it. This process is repeated in a limit cycle (Muller et al. 2003; Karpen 
et al. 2003) with a period of ~ 3 hrs. In this cycle a pair of consecutive condensations forms, 
moves along the tube, and finally falls to the chromosphere. The difference between tubes 
1 and 2 in terms of the plasma evolution is dictated by the heating asymmetry relative to 
the flux tube geometry (see, e.g., Karpen et al. 2003, 2005). In both cases the heating is 
stronger in the righthand footpoints of Figure 2. In tube 1 the steep slope between cl = 40 
Mm and 110 Mm keeps the condensation in the dip. However, in tube 2 the slope is not 
steep enough to retain the condensation, so it falls to the lefthand chromosphere. 

The evolution of the mass in shallow-dip tubes, simulations 3 and 4, is shown in Figure 
4b. In these two simulation cases there are no stable condensations as in simulation 1; only 
small and short-lived condensations are produced cyclically. In simulation 3, the sequence of 
condensing and falling to the chromosphere is repeated every 2.5 hrs, and each condensation 
lives for approximately 0.5 hrs. However, in simulation 4 the plasma condenses in more 
massive condensations roughly every 4 hrs, and each has a lifetime of 3 hrs. Tube 3 has a 
length of 176 Mm, and tube 4 has a length of 240 Mm. This indicates that the longer the 
loop, the larger the period of the limit cycle. 

Similarly, the evolution of the cool mass of the two loop-like tubes is plotted in Figure 
4c. Small condensations form in both loops of the overarching arcade. These two tubes (5 
and 6) have lengths of 121 Mm and 479 Mm respectively, with tube 6 being the longest tube 
considered in this work. In simulation 5, we see that each condensation forms quickly, lives 
for ~ 15 min, then falls rapidly to the chromosphere. This process is repeated periodically 
every 3 hrs. The temporal evolution of the cool mass in loop 6 is more complex, with several 
peaks reflecting a cluster of short-lived condensations that forms and disappears within an 
interval of ~ 6 hours. Only one entire condensation cycle occurs within the simulation 
time, but another condensation appears at the end of the run. Thus, the period of the 
evaporation/condensation cycle is around 18 hrs. 


3.1. Condensation masses and lengths 

Figure 5 shows the maximum mass of the condensation in each flux tube. This quantity 
is useful for determining the ability of the tube to produce and accumulate cool mass accord- 
ing to its geometry. When the condensation mass grows monotonically, as in simulation 1, 
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Fig. 5. — Plot of the maximum mass of the condensation as a function of tube length. 
Diamonds symbols correspond to condensations of the population one with masses larger 
than the separator mass (dotted line). Asterisks corresponds to condensations of population 
two with masses smaller than the separator. The cases of the representative tubes are 
surrounded with black circles (see text). 

this maximum mass coincides with the value at the end of the simulation. When N consec- 
utive evaporation/condensation cycles are produced, as in simulations 2 to 6, we average the 
maximum mass at each cycle as M max = where m™ ax is the maximum mass at 

i-th cycle. From this figure we see that the condensations can be divided mainly into two 
separate groups: condensations of relatively large masses (diamonds), and condensations of 
relatively small masses (asterisks), ft is important to note that all 125 simulated flux tubes 
studied here produce condensations. In the first group, condensations form and then grow 
monotonically until the end of the simulation, as for simulation 1 (see Fig. 2). In the second 
group, each condensation forms, moves along the tube, and finally falls to the chromosphere, 
as in tubes 2 to 6. We have found that approximately 7 x 10” 4 g cm -2 is the separator be- 
tween the two condensation populations (dotted line). Condensations just above this limit 
initially have cycles of evaporation/condensation and fall to the chromosphere, but eventu- 
ally the condensation mass grows monotonically. Condensations just below the limit start to 
grow monotonically but at some time the entire condensation falls to the chromosphere, and 
several cycles are produced in the simulation. All flux tubes that produce condensations of 
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the first group are tubes with deep dips, but not all tubes with deep dips produce conden- 
sations of the first group, as discussed above for tube 2. Small condensations of the second 
group are produced in all flux tube geometries, as the representative examples 2 to 6 indicate. 
There is no clear dependence of the condensation mass on the tube length. Condensations 
of the first group are produced in a range of tube lengths between 160 and 340 Mm. This 
reflects the range of lengths of the deeply dipped held lines (see Fig. 3). In Figure 5 we see 
a few examples of the second group of condensations for tube lengths larger than 350 Mm. 
These condensations are generated in the overlying arcade, and as mentioned in §2 we have 
only a small sample of those held lines. However, for our small sampling of the overlying 
arcade, the maximum condensation mass does not depend on tube length. In summary, all 
three kinds of tube geometries produce condensations of the second population, with small 
masses, but a subset of deeply dipped hux tubes yields the massive condensations of the hrst 
population. 

In Figure 6, the mass of the hot plasma and the mass of the condensations of all 
tubes are plotted. The mass of the hot plasma (dotted line) grows steadily from a small 
value (25 x 10 -4 g cm -2 ) given by the initial uniformly heated equilibrium. In this interval 
the initially underdense tubes are being filled with hot plasma as a consequence of the 
chromospheric evaporation at the tube feet. The growth of the hot plasma slows down at 
about t = 2.0 hrs, and finally stops at t = 2.5 hrs (see inset plot of Fig. 6). In the interval 
2. 0-2. 5 hrs, the chromospheric evaporation has increased the coronal density of most of the 
flux tubes and also their radiative losses, so the tubes start to cool down. At t = 2.0 hrs, the 
hrst condensations appear and the mass of the cool plasma starts to grow. It is interesting 
to note that in the interval 2. 5-4.0 hrs the coronal mass is reduced by 30% of its maximum 
value. This partial evacuation of the corona occurs because the cool condensations accrete 
mass faster than the process of evaporation from the footpoints, due to a pressure deficit 
that pulls additional material into the condensations from the corona (see Karpen et al. 
2003). The growth rate of the condensation mass is higher during this interval than in the 
rest of the simulation. After this time ( t = 4.0 hrs) the hot mass reaches a more or less 
constant value of 158 x 10 -4 g cm -2 , on average, for the rest of the simulation. Similarly, 
the condensation mass reaches a steady growth rate of 39 x 10~ 4 g cm -2 hr -1 . In this steady 
state, the rate at which mass condenses is balanced on average by the rate of evaporation 
from the footpoints of all tubes, and the mass of the coronal plasma remains more or less 
constant. 

We have separated the contributions of the populations one (dashed line) and two (dot- 
dashed line) to the total cool mass. Figure 6 shows that there is only a slight difference 
between the condensation mass of the entire system and the mass in tubes of population one. 
This indicates that the main contribution to prominence mass comes from the population 



13 



Fig. 6. — Temporal evolution of the mass of the plasma in the filament channel. The mass of 
the condensations (solid line) is separated according to the contributions of population one 
(dashed line) and population two (dot-dashed line) condensations. The mass of the coronal 
plasma (dotted line) grows from an initial equilibrium value beginning at the localized heating 
onset (t — 0 hrs). The inset plot shows details of the early evolution. 

one tubes. Lin et al. (2005) estimated the observed width of the prominence threads to 
be 200 km. Assuming this value for the average diameter at the deep-clip tube centers, we 
can estimate the total mass of our model prominence to be 1 x 10 14 g. This agrees with the 
observed values ranging from 1 x 10 14 g to 2 x 10 15 g (see Labrosse et al. 2010, for an extended 
review and references therein). The mass estimate depends strongly on the volumetric filling 
factor of the cool plasma within the surrounding hot corona, which can be estimated by 
dividing the coronal volume occupied by cool threads by the total volume in which those 
threads are sparsely distributed. This total volume was determined as described in §4 and 
is ~ 23 x 10 3 Mm 3 . The 47 population one threads collectively occupy ~ 54 Mm 3 assuming 
cylindrical threads of 200 km in diameter and the average thread length at the end of the run 
derived from Figure 8. Therefore the filling factor of our model is ~ 0.002, which matches 
the lower limit of observed values (Labrosse et al. 2010). To obtain a larger filling factor, as 
reported from some observations, either the threads are thicker than assumed (which would 
contradict the highest resolution data) or the prominence contains appreciably more threads 
than were selected for our modeling investigation. 
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Fig. 7. — The maximum mass of every individual condensation (see Fig. 5 for definition) 
displayed versus its length. 

In Figure 7, we have plotted the maximum condensation mass as a function of its 
maximum length. We see that condensations of the first population are long, whereas con- 
densations of the second population are mainly short. This result is expected because when 
a condensation grows by accreting mass, it also grows in size. For the first population there 
is no clear relation between the mass and length of the condensations. Lengths are in the 
range 26.4 to 56.6 Mm, with an average value of 37 Mm. These large condensations form 
threads similar to those reported by Lin et al. (2005); the computed lengths are in agreement 
with the observed lengths. We hereafter call the condensations of population one “threads” . 
Condensations of population two are smaller, with an average maximum length of 3.3 Mm. 
Hereafter, these small features are called “blobs” . There is observational evidence for these 
features in the corona, both in prominences and in coronal rain (e.g., Schrijver 2001; Muller 
et al. 2003, 2004; De Groof et al. 2004; de Groof et al. 2005). 

In Figure 8, the averaged lengths of the two populations of condensations are plotted. 
The size of the threads (solid line) grows rapidly between t = 2.0 and t = 3.5 hrs, when 
the average size reaches 8 Mm. This period of time coincides with the partial evacuation of 
the coronal mass discussed in §3.1. Note also the sudden increase of the mean size between 
t = 3.0 hrs and t = 3.3 hrs; in this 18 min interval the mean length grows from 2 to 8 Mm. 
Later on, the length remains more or less constant at about 8 Mm, until i « 5 hrs when 
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the growth resumes up to the end of the simulation. In contrast to the total mass (Fig. 
6), the length does not grow linearly, and its growth rate decreases with time. Two factors 
contribute to the diminution of the growth rate, and both are related to the geometry of 
the deep-dip flux tubes. First, when the condensation grows the newly accreted material 
occupies a region of the tube with higher gravity due to the curvature, which compresses the 
condensation more dne to the small gravitational scale height. Second, the newly accreted 
material fills regions of the tube that are wider dne to the area expansion. In contrast to the 
threads, blobs rapidly reach a constant length of about 1 Mm. These small condensations 
can form anywhere in the tube and are short-lived (see Fig. 4). This continual process of 
condensation and falling to the chromosphere keeps the temporally averaged length more or 
less constant. 



Fig. 8. — The averaged length of the condensations in each population as a function of time. 


3.2. Condensation velocities 

After formation a condensation moves along its held line. Its motion is determined 
by the local gravity along the tube and the localized-heating imbalance between the tube 
footpoints. Condensations are extended regions of cool and dense material in which each 
point moves with its own velocity. These motions can be decomposed into a bulk motion of 
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Fig. 9. — Absolute values of the maximum bulk speeds of the condensations of both popu- 
lations as a function of the angle between the bulk velocity vector and the vertical 2 -axis, 6. 
The angle ranges from 9 — it/ 2 (horizontal motion) to 9 = tt (vertical downward motion). 

the condensation plus local expansions or contractions. Here we investigate the dynamics of 
these individual condensations considering the bulk motion as the mass-averaged velocity in 
the condensation, i.e., the motion of the center of mass. In Figure 9, the absolute value of 
the maximum bulk velocity of the condensations is plotted. This magnitude is shown as a 
function of the angle between the bulk velocity vector and the vertical direction. There are 
two kinds of motions associated with the two condensation populations. Thread maximum 
velocities are concentrated in an angle range 9 = 7t/ 2-2.0 that is almost horizontal. Figure 9 
reflects the fact that motion of the threads is in the concave part of the tube, and the motion 
in the dip is mainly horizontal with a small vertical projection, resembling a swinging motion. 
The direction of the velocities is not included in the figure, but in the movies associated 
with Figures 10 and 11 we see horizontal counter-streaming in adjacent threads. These 
condensations form in the dipped part of the tube, and fall to the bottom of the dip with 
a damped oscillatory motion (Antiochos et al. 2000; Karpen et al. 2003). The maximum 
velocity is reached beyond the bottom of the dip; this is why the angles are not exactly 7 t/ 2 
in Figure 9. In contrast, blob velocities are in a range 9 = n/2 — 2.9, indicating that blobs 
move from near horizontal direction, 9 = n/2, to near vertical downward direction, 9 = it. 
This result is expected because blobs are produced on flux tubes of all three kinds studied 
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in §2, with geometries ranging from mostly horizontal tubes to mostly vertical tubes. The 
motion of the blobs on tubes with either deep dips or shallow dips (simulations 2 to 4) is 
relatively horizontal, while that of blobs formed in the overlying arcade (simulations 5 and 
6) is nearly vertical. Most of the blobs slow down as they approach the tube feet, and in 
some cases several bounces occur before they fall finally to the chromosphere (see movies 
associated with the Ha proxy of Figs. 10 and 11). In general the maximum velocities are not 
reached at the tube feet. This behavior has been observed by Schrijver (2001) and studied 
theoretically by Muller et al. (2003). Threads have a maximum velocity of 59.9 km s^ 1 
and a minimum of 6.3 km s _1 , with a mean value of 28.8 km sW Blobs have a maximum 
velocity of 83 km s^ 1 and a minimum of 14.6 km s -1 , with a mean of 45.9 km s -1 . There is 
no clear dependence of the bulk velocity on the angle in the thread population. However, 
the blobs show a small tendency for the maximum bulk velocity to increase with the angle. 
This indicates that the downward velocities are larger in vertical tubes, but these velocities 
are far from the free-fall speeds. As we see in Figures 10 and 11, the blobs form mainly 
below 30 Mm in height. At this height the free-fall velocity is approximately 130 km s^ 1 
(Schrijver 2001; Mackay & Galsgaard 2001), indicating that the condensations are sustained 
by pressure forces impeding the free fall. 


4. Observational diagnostics 

As in our earlier works on thermal nonequilibrium, we illustrate the dynamic behavior 
of the entire prominence model by computing synthetic images in selected representative 
EUV passbands, using the temperature response function of the Atmospheric Imaging As- 
sembly (AIA) instrument aboard the Solar Dynamics Observatory (SDO) satellite. Unlike 
our prior studies of individual tubes seen from a single point of view, these images reflect the 
3D structure of the prominence model with orientation and line-of-sight integration effects 
included. We show the AIA band passes 17lA, 21lA, and Ha proxy. Several spectral lines 
contribute to each EUV passband, associated with different ions and ionization degrees. The 
peak temperature of the 17lA passband is log T = 5.8, and the main contribution comes 
from Fe IX. The 21lA band has a peak temperature of log T = 6.3, and the main contribu- 
tion is from Fe XIV. We have assumed that the emission in the EUV channels is optically 
thin and computed the images using the temperature response function of the SDO/AIA 
instrument. The Ha emission is computed as in our earlier papers (see, e.g., Karpen et al. 
2001), setting the intensity constant whenever T < 35,000 K and zero at higher temper- 
atures. We assume that the flux tubes are unresolved by the telescope, and use the AIA 
instruments point spread function half-width of 0.6”. 



We have constructed off-limb images in two perpendicular orientations of the system 
relative to the line of sight (LOS). In the end-view orientation the LOS is parallel to the PIL 
or filament-channel axis. In the side view the relative orientation of the LOS is perpendicular 
to the PIL. Figures 10 and 11 show the temporal evolution in all of the considered channels 
in the end and side views, respectively. In both views different types of structures are visible 
in different emission lines. In the Ho; channel we see the temporal evolution of the cool 
condensations, and distinguish between the two condensation populations. As we described 
in §3.1 the threads are the largest cool condensations, containing most of the prominence 
mass. From the Ha images we identify two separated prominence segments in the sheared 
double arcade. Both have a trapezoidal shape with a base of approximately 20 Mm width 
(end view), 115 Mm length (side view), and 20 Mm height. Threads form in an off-center 
position in the dipped portion of the flux tubes, and oscillate around the bottom of the 
dips. In Figure 11a the recently formed threads are shown. The average length of the 
threads is approximately 10 Mm (see Fig. 8), and they move mainly horizontally with a 
large amplitude. In the movies accompanying the figures we clearly see the opposite motions 
on adjacent threads. Antiochos et al. (1999) proposed this mechanism of condensation and 
oscillation in the dips of the flux tubes to explain the counterstreaming motions reported by 
Zirker et al. (1998). Here we clearly see how a bundle of threads manifests counterstreaming 
motions. These motions are also present in the end-view orientation (Fig. 10b). However, 
the direction of the motion is slightly different with respect to the LOS and the projected 
amplitude of motion is small. In Figure lib the threads are larger than in the previous 
frame 11a, with an average length of 35 Mm. We clearly see the concave- up shape of the 
threads in Ha (Figures 10b and lib) lying in the bottom of the dips. The amplitude of the 
oscillations is now much smaller than in the initial stages of the oscillation, because of the 
strong damping reported by Antiochos et al. (2000). 

Condensations of the second population (blobs) are small features seen in the Ha syn- 
thetic images in all regions of the filament channel and the overlying arcade. These ubiquitous 
features form mainly at heights below 50 Mm, although some blobs form at greater heights 
up to 80 Mm. Blobs forming in the overlying arcade rapidly fall to the chromosphere with 
an accelerated motion following the magnetic held. As we pointed out in §3.2, these motions 
are not free-fall movements. When blobs approach the tube footpoint, their movement decel- 
erates before they enter into the chromosphere. Schrijver (2001) reported similar motions in 
the coronal rain observed by TRACE, and suggested that the deceleration is the consequence 
of the pressure gradients of the underlying plasma close to the tube footpoints. Muller et al. 
(2003) simulated this phenomenon and confirmed that blobs increase the pressure of the 
plasma below as they approach the chromosphere. 

Figures 10c and lOd show synthetic images in the 17lA passband in the end-view orien- 
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Fig. 10. — Time evolution of the whole system in end-view with different filters. In (a) 
and (b) the Ha proxy is shown at t — 5.56 hrs and t = 23.81 hrs, respectively. Similarly, 
(c) and (d) show these two stages of the evolution in the 17lA passband of the AIA/SDO 
instrument. In (e) and (f) the temporal evolution in the 21lA emission intensity is plotted. 
(An animation of this figure is available in the online journal.) 











tation. In this channel the emission comes from relatively warm plasma at log T = 5.8. In 
these figures, we clearly see strong emission from the core of the filament channel, where the 
prominence resides. The source of this emission is the warm plasma at the prominence-corona 
transition region (PCTR) at both ends of the individual threads (see the corresponding Ha 
images). The emission in this wavelength is strong in the core because the threads are con- 
centrated in this region and the bright features are more or less aligned with the LOS. In 
the side views (Figs. 11a and lib) the emission from the filament core is not so strong, be- 
cause the bright features are less concentrated and not aligned with the LOS. However, from 
this point of view we clearly see the bright structures resembling two opposing tadpoles (or 
comets with tails) on both ends of the threads. Between each pair of tadpoles is an emission 
gap associated with the cool material of the condensation. These tadpole features are the 
consequence of the temperature and density distributions at the ends of the condensations. 
The temperature rapidly rises from several thousand Kelvin near the condensation to sev- 
eral million degrees in the hot parts of the tube, while the density inversely declines. In an 
intermediate region between the condensation and the corona, the temperature and density 
conditions are appropriate to produce the bright emission. The gradients of the tempera- 
ture and density steepen close to the condensation and flatten farther away, explaining the 
tadpole shape of the emission. These bright features move together with the condensations, 
and as we show in the accompanying mpg movies, counterstreaming motions then are also 
evident in the 17lA passband. The end-view images reveal that the emission from the outer 
region of the prominence is fainter than from the filament core. In this region, the bright 
features are associated with condensations of population two blobs. Since blobs are small 
condensations, the emission gap between their pairs of tadpoles is smaller than in case of the 
threads. The emission depends on the motion of the blobs: the leading tadpole is brighter 
than the trailing tadpole, relative to the direction of the blob motion, and in some cases the 
trailing tadpole is very faint. The animations also reveal long and transient brightenings, 
which in some cases illuminate most of the tube. These very long, bright features are the 
evaporation/condensation cycles. Before a condensation is produced, a large portion of the 
material of the flux tube cools down and the density increases. When the temperature of this 
denser plasma is approximately the peak temperature of the filter, a strong brightening is 
produced in that passband. The plasma continues to cool down, and finally the condensation 
is produced. At this moment the long, bright feature disappears and the tadpoles appear in 
17lA. All of the flux tubes in the filament channel produce such long features. However, the 
tubes with the population-one condensations produce this transient brightening only in the 
initial stages of the simulation, just before thread formation (around t — 2 hrs). Thus, we 
expect large transient brightenings to occur in the 17lA passband, both in the core of the 
FC and in the overlying arcade. 
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Figures lOe and lOf illustrate the temporal evolution of the system in the end-view 
orientation for the hottest channel, 21 lA. In these images we see that hot plasma fills most 
of the FC structure, without the small features that dominate the 171 A images because 
the 21lA channel has a wide temperature response function (see, e.g., Lernen et al. 2011). 
The hot material is evaporated from the chromospheric plasma by the footpoint heating. In 
both images, the core of the structure is bright, resembling the “chewy nougats” observed 
by Yohkoh (Hudson et al. 1999; Hudson & Schwenn 2000). In contrast with the 17lA 
images, however a diffuse halo surrounds the prominence in 21 lA emission reaching a height 
of 40 Mm, resembling the “hot prominence shrouds” reported by Habbal et al. (2010). In 
the side views (Figs, lie and Ilf) the emission is more uniform, with no distinguishable 
bright structure in the core. This demonstrates that the intensity distribution depends on 
the relative orientation of the FC to the LOS. In particular, the bright core appears when 
the PIL is nearly aligned with the LOS. In the end view, the flux tubes in the core are closely 
aligned with the LOS, producing the strong core emission. In the side views (Figs, lie and 
Ilf), the 21lA emission is more uniform and the core is darker because the cool plasma of the 
prominence emits little in this bandpass, an effect also known as emissivity blocking (Anzer 
et al. 2007). In addition, the emission in this channel is more or less constant throughout 
the simulation time (see the accompanying movies) in both orientations. This indicates less 
dynamism of the hotter plasma than the cooler plasma. 


4.1. Differential emission measure 

In Figure 12, the differential emission measure (DEM) of our prominence model and 
overlying arcade is shown, taking into account all the flux tubes of the filament channel 
structure. This is averaged in space and time for each flux tube, and then combined over 
the whole volume. To avoid the contribution of the initial transient phase, we averaged the 
DEM from 5 hrs after the onset of the localized heating to the end of the simulation. The 
DEM is computed as in a previous paper (Karpen & Antiochos 2008). The result of this 
computation is normalized by a prescribed column depth, which is chosen to ensure a good 
fit with the DEM values derived from prominence observations by SOHO SUMER (Parenti 
& Vial 2007). The resulting column depth is h ~ 10 km, two orders of magnitude below our 
previous results ( h ~ 1000 km) because the present work includes approximately 100 times as 
many tubes. As we see in the figure, our model reproduces very well the shape of the observed 
quiet- prominence DEM between log T = 4.6 and 5.7. In addition, the location of the relative 
minimum of our profile is log T = 5.1, consistent with the observational value. As argued 
in Karpen & Antiochos (2008), the additional emission at higher temperatures (log 7 « 6) 
can be attributed to the unsubtracted background and foreground coronal contribution to 
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the observed DEM of Parenti & Vial (2007). Moreover, in our filament channel model we 
have considered only a small fraction of the overlying arcade loops that exist in a realistic 
structure. These loops of the overlying arcade contribute mainly to the warm and hot part 
of the DEM function because the condensations formed in those loops are small. Therefore, 
increasing the number of overlying arcade loops should yield better agreement between the 
computed DEM and the observed curve at temperatures higher than log T = 5.7. We plan 
to extend our model to include more of the overlying arcade. 


5. Discussion and conclusions 

We have developed a comprehensive three-dimensional model of the plasma dynamics 
and energetics of a prominence and its overlying arcade. We have found that condensations 
are formed in all of the flux tubes studied. There are two populations of condensations 
characterized by their cool mass. If the cool mass is below 7 x 10” 4 g cm -2 , repeated 
cycles of condensations forming and falling to the chromosphere are produced. Above this 
limit, condensations grow at a constant rate, reaching large values for the mass and length. 
We have called the small condensations “blobs” and the large ones “threads”. Blobs are 
small, with an average length of about 1 Mm, while threads are larger condensations with 
an average length of 37 Mm. The dominant contribution to the prominence mass comes 
from the threads. This mass grows at a constant rate of 39 x 10” 4 g cm -2 hr” 1 . Assuming 
that the threads have a width of 200 km, the total mass of the prominence is 1 x 10 14 g at 
the end of the simulation, in agreement with observed values. When threads are formed, 
they oscillate around an equilibrium position; this motion is mainly horizontal (parallel to 
the solar surface), with a mean velocity for the entire set of threads of 29 km s” 1 and a 
maximum value of 60 km s” 1 . Blobs move in all directions from horizontal to vertical, with 
a mean velocity of 46 km s” 1 and a maximum of 83 km s” 1 . Blobs of the overlying arcade 
travel downward, and are the coronal rain observed in coronal loops (Schrijver 2001; De 
Groof et al. 2004; Mlillcr et ah 2003, 2004, 2005). 

With our three-dimensional model we have constructed synthetic images including the 
orientation of the system and LOS integration effects. The synthetic Ha images reveal 
two trapezoidal-shaped prominence segments, each centered at a magnetic dipole and made 
up of threads. These images clearly show the threads growing and the counterstreaming 
motion occurring on adjacent threads. We have found that blobs are ubiquitous features 
in these images, moving along the held lines and falling to the chromosphere. Most of the 
blobs decelerate before entering the chromosphere. Blobs of the overlying arcade are the 
coronal rain, with primarily downward motions. In the sheared flux tubes associated with 



23 


the prominence, blobs are also produced, but their motions are more horizontal than the 
coronal rain. We have also computed synthetic images in two EUV channels, 17lA and 
21lA. In 17lA, we predict bright emission from the core of the FC where the prominence is 
located. The emission in that region comes mainly from PCTR warm plasma. The region 
surrounding the core is darker, but transient moving features appear whose emission comes 
from the warm plasma between the cool blobs and the hot corona. In addition, there are long 
bright transient features occupying most of a tube, forming just before a blob is produced. 
In 211 A, strong emission comes from the FC core in the end-view orientation, and a diffuse 
halo surrounding the prominence. This indicates that the cool material of the prominence is 
surrounded by hot material resembling the “chewy nougats” seen by Yohkoh (Hudson et al. 
1999) and the observed “hot prominence shroud” from recent eclipses (Habbal et al. 2010). 
We have computed the DEM and found very good agreement with the observational DEM 
of a quiescent prominence (Parenti & Vial 2007). The computed DEM curve reproduces the 
observed shape between temperatures log T = 4.6 and 5.7, and the observed position of the 
minimum at log T = 5.1. 

The ability of our model to reproduce many observed prominence characteristics lends 
strong support to the thermal non-equilibrium mechanism for producing the mass of the 
prominence, as well as the plasma dynamics. The agreement between model and observations 
also suggests that the fundamental structure of the filament channel is well described by the 
sheared arcade. Currently, we are extending the realism of our model in three directions. 
First, we are greatly increasing the number of tubes that sample the overlying arcade rooted 
outside the filament channel. Second, in this work have assumed that the footpoint heating 
is steady in time and is the same on all flux tubes. In an upcoming work we will include the 
dependence of the heating function on the key physical parameters, according to the scaling 
laws summarized in Mandrini et al. (2000), as well as allow for temporally varying heating. 
Finally, the emission in EUV coronal lines is significantly affected by Lyman absorption by 
neutral H and He and by He + (e.g., Kucera et al. 1998); these important effects will be 
included in future visualizations of our simulations. The large amplitude oscillations of the 
condensations in our prominence model will be addressed in an upcoming publication. 
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Fig. 11. — Same as Fig. 10 showing the temporal evolution in emission intensities in the 
three filters, but from the side-view orientation. (An animation of this figure is available in 
the online journal.) 










29 



Fig. 12. — Time-averaged differential emission measure vs. temperature for the whole system 
of flux tubes of the filament channel and the corresponding overlying arcade (solid line). The 
observed best fit obtained by Parenti & Vial (2007) for prominence section A_1 (data courtesy 
of Parenti 2007 [private communication]) also is shown (dashed line). 




